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I. ION OPTICS 


A major effort has been undertaken in ion op- 
tics, with specific emphasis on theoretical calcu- 
lations and the construction of an inhomogeneous 
field instrument. A 255° magnet of large dimensions 
was conceived in order to provide a laboratory in- 
strument of great versatility. 

As a mass spectrometer, this analyzer which is 
being developed should have the maximum sensitivity 
of any single stage analyzer. The applications of 
a single stage, high- transmission, double focusing, 

large radius analyzing system are virtually limit- 

12 3 

less. ^ ’ Lunar, meteoric and geological samples 
as well as medical tracer specimens could also be 
analyzed for isotopic content or low level alpha 
activity, (see Table A). Analyses, similar to the 
one reported by Mewherter et al^-V could also be 
made for plutonium or any other man-made isotope 
which may occur in the environment. Because of the 
high transmission characteristic, the thermal ion- 
ization efficiencies of various source materials 
can be measured in addition to providing experi- 
mental data for determining work functions. Both 
negative and positive ionization efficiencies could 
be measured by reversing the analyzing magnetic 
field. In the study of crystal structure, surface 
phenomena, and ion transmission properties, the 
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ALPHA PARTICLE ACTIVITIES OF THE TRANSURANIC ELEMENTS 
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Specific activities times % of emission of the stated alpha line Table 


magnet could provide a basic analyzing system. The 
large radius will allow a particle energy range of 
10 percent while a 4 inch gap width permits installa- 
tion of crystalline target material within the mag- 
netic field. High energy a-particle backscattering 
experiments or improved ion-transmission systems 
such as the one described by Rasekhi^ can be conducted 
with this device. Field emission ion sources or even 
an electric lens for improved resolution may be in- 
stalled within the magnetic field. Since this device 
can accommodate high energy ions, sputtering sources 
with high potential acceleration anodes may be 
utilized. The large energy spread of such ions great- 
ly affects the resolving power of an anlyzer. How- 
ever, by accelerating the sputtered ions before 
analysis, their relative energy dispersion is reduced. 
Low intensity processes such as photo-ionization would 
be well suited for investigation by this high-trans- 
mission system. Or conversely the analyzer could 
provide an intense ion beam of various energies for the 
study of ion to photon conversion properties. The 
large radial and axial extent of the magnetic field may 
make it possible to develop s point to point ion- 
imaging system in order to investigate spatial surface 
emission characteristics of a sample. 

The energy of a-particles, like any other charged- 
particle radiation, can be measured by the number of 
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ion pairs created by it upon passing through a me- 
dium. Some individual particle detectors are: 
ionization and cloud chambers, proportional, spark, 
and scintillation counters, nuclear emulsions, and 
semiconductor devices. Only the last three are 
easily adapted to charged-particle spectrometry. 
Nuclear emulsions developed for high energy experi- 
ments and spectrographic detectors are highly 
efficient and extremely sensitive. Such plates are 
useful for recording rare events due to their in- 
sensitivity to background. However, it is impos- 
sible to separate events in time and the resultant 
energy resolution is low ~5%. Scintillation^ coun- 
ters are simple, reliable,' and have fast response 
times as well as a large working area and insensi- 
tivity to intense beta and gamma background radia- 
tion. Again, these devices exhibit a low energy 
resolution ~57o. Surface barrier detectors are per- 
haps the most versatile radiation detectors commonly 

used. Nevertheless, their energy resolution and 

/ 

background level are inferior to magnetic spectro- 
graphs. Combining a high-transmission magnetic 
analyzer with a solid-state detection system would 
quite probably result in an optimum analyzing system, 
superior to either device operating independently. 
Range and ionization fluctuations inherent in the de- 
tector would be minimal. Such a system would be well 
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suited to detect low level alpha-radiation emitted in 
the presence of competing strong reactions which pro- 
duce high background radiation. 

Luminosity and resolution are the primary requi- 
sites for selecting a magnetic optical system for the 
charged-particle investigations described previously. 
High luminosity is characterized by a large acceptance 
angle and the 255*^ system has a greater accep- 

tance angle (~0.37o of 4tt steradians) than an equiva- 
lent 180° device. One characteristic of 180° inhomo- 
geneous magnets is that their ion source and focal 
point are in the fringing fields of the system. With 
a 255° device, it is possible to locate both the 
source and detector within the main field of the magnet 
(Fig. 1). Several 180° instruments have been con- 
structed for investigating a-particle radiation, but 
at present there are no 255° devices under construc- 
tion in this country for use in the research activities 
previously mentioned. 

Dzhelepov et al.,^ have constructed a 255° device 
promarily for the investigation of a-particles. The 
estimated resolution is ~0.10% with a transmission 
factor approximately 0.3% of 4 tt steradians. This 
effective solid angle is at least an order of magnitude 
greater than other types of spectrometers. Photo- 
graphic emulsions are used in the detection system even 
though solid state detectors would allow direct, real- 
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time readout capability. Shaping of the pole pieces 
was accomplished by field expansion techniques and 
experimentally determined mechanical shims. Experi- 
ments conducted with this device indicate that it 
operates within the theoretically predicted resolu- 
tion and transmission parameters. 

Several 180° inhomogeneous field spectrometers 

8 9 

such as those by Holtebekk and by Judd have been 
constructed for charged-particle investigations. 

_3 

Acceptance solid angles are quoted at about 1.2’ 10 
of 4n steradians with a resolving power of approxi- 
mately 1000. Field expansion techniques were utilized 
in the magnetic field design and 3rd order coeffi- 
cients were calculated in order to reduce aberrations. 
Preliminary measurements on the Judd instrument at 
Berkeley indicate a resolution of 0.17o with a solid 
angle of about 7*10”^ of 4 tt steradians. 

Specifications of Analyzing Magnet 

The RPI Analyzing magnet has been constructed 
to the following parameters: 

1.0 Mean radius of 27.8 inches 

2.0 Maximum extent of pole pieces is 270° 

3.0 Total beam deflection angle is 255° (V^2n) 

4.0 Field Index of inhomogeneity is 0.5 

5.0 Gap width is 4.0 inches at the mean radius 

6.0 Pole pieces extend from 16 to 44 inches 

7.0 Maximum induction in the gap is 5500 gauss 

at a radius of 27.8 inches 
7 



8.0 Maximum magnetic induction in the pole 
pieces is 11,200 gauss 

9.0 Object and image points are at the mean 
radius within the field 

10.0 Design beam acceptance angle is + 15° 
in the horizontal plane and + 2° in the 
vertical plane 

11.0 Total acceptance of 0.3% of 4n steradians 

12.0 Resolving power of approximately 1000, 
with maximum beam acceptance 

13.0 Magnet weight is 70 tons. 

Excitation is provided by four exciting coils conform- 
ing to the following properties: 

1.0 Two copper conductors 0.530” x 0.100” 
in parallel 

2.0 Pyre ML film insulation on all conduc- 
tors 

3.0 Coils constructed on 88 inch diameter 
aluminum flanges 

4.0 Water cooling around the edge of the 
coil with a total flow of 1.73 gal. /min. 
through 0.25 inch tubes resulting in a 
20 degree Fahrenheit temperature rise 

5.0 Coil resistance (two wires in parallel) 
is 2.53 ohms 

6.0 Weight of copper is 1900 lbs.; weight 
of aluminum flange is 138 lbs. 
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7.0 Coil current is 23.7 amperes at 60 
volts . 

8.0 Magnetic excitation is 13,050 ampere 
turns/coil 

9.0 Power dissipation for each coil is 
1,422 watts 

Construction of the iron yoke is complete and is pic- 
tured in Fig. 2. The assembly consists of three 
legs and six yoke pieces as depicted. The coil assem- 
blies are shown in Figure 3 and a constructed working 
model in Fig. 4. Figure 5 is a photograph of a pole 
piece in. the final machining stage. Power, vacuum, and 

instrumentation is as yet incomplete. 

Theoretical Considerations of Sector Magnetic Fields 
The RPI analyzing magnet is a 255° device exhib- 
iting azimuthal symmetry about the z axis and mirror 
symmetry across the magnetic midplane. Figure 6 de- 
picts the general pole piece shape, utilizing an in- 
homogeneous field concept. The magnetic field B is 
a decreasing function of the path radius r. The 
forces exerted on a particle traversing this field 
with velocity v are indicated in the lower quadrant. 

The radial and axial forces tend to restore the 
particle to its initial values of r and z. This prop- 
erty is characteristic of inhomogeneous field systems 
and double focusing (radial as well as axial) is 
possible. The 255° analyzing lens has this property 
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Exciting Coil 
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Figure 4 Model of 255 Inhomogeneous Magnetic Alpha Particle and Mass 

Spectrometer Assembly 
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Figure 6 
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and with reduced transmission, there is a point to 
point object-image correspondence. Figure 7 depicts 
the axial focusing properties of a 255° analyzer, 
while Fig. 1 indicates the radial properties. The 
focal properties of this system are completely deter- 
mined by the midplane induction field expansion, 
which is a Neumann boundary condition. Unfortunately, 
it is not possible to choose these coefficients such 
that the optimum radial and axial focusing properties 
are achieved simultaneously.^^ Since the most impor- 
tant characteristic of an ion analyzer is resolution 
as opposed to minimiam focal area, this system will 
be designed to minimize radial image divergence. 

In a cylindrical coordinate system, the magnetic 
field at the analyzer midplane can be represented by 
the function: 

where 

Bq = the magnetic induction at r = r^ , z = 0 

r = the mean radius 
o 

n = the inhomogeneous field index 
In general, the field index will have a value between 
0.0 and 1.0 for optimum radial focusing. The magnetic 
field, B, can be derived either from a vector field A 
or a Scalar magnetic potential Vm. For the development 
and solution of the general equations of motion, the 
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Figure 7 
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vector potential A is appropriate, however, the 
scalar potential Vm is necessary for determining 
isomagnetic potential contours. Both approaches 
are straight forward and only algebraic relation- 
ships between field expansion coefficients are 
required to transfer from the vector field deriva- 

— k 

tion of B to the magnetic potential description of 

B. 

Derivation of the Midplane Magnetic Field from a 
Vector Potential and a Potential Function 

The magnetic field induction across the midplane 
can be expanded in a Taylor series about r^ as follows: 

= BQ[l+bj^P+b2P^+b^p^+b^p^+b^p^+ ...] 

where 



bj^, b2 = the field shaping constants. 

The first shaping constant is related to the basic field 
inhoraogenity value by 



b2, b^, b^ are higher order aberration parameters 

which are varied to obtain the desired focal properties 
These five coefficients uniquely determine the magnetic 
field in three dimensions, through fifth order. 

•A 

B may be derived from a vector potential by 
B = V X A . 
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For a cylindrical coordinate system, 

R = _Ji + -L A6 

r^ 5p r^ (1+p) 

and since the magnetic field is azimuthally invariant, 

Aq = B^r^(p+a2p^+a2p^+a^p^+a^p^+a^p^. . ) 
we have 

[Bz]z=o = B^[l+(l+2a2)P + (-l+a2+3a3)P^ 

+ (l-a2+23+4a^)p^ + (-l+a2-a3+a^+5a3)p^ 

+ (l-a2+a3~a^+a3+6a^) p^+. . ] . 

Comparing coefficients of equal powers of p yields 
= 2a2 + 1 
b^ 3a^ ^ ^ 

b^ “ ^3 "* ^2 ^ ^ 

^4 = 5a3 + a^ - a3 + a2 - 1 

= 6a^ + a^ - a4 + a3 - a2 + 1 • 

is determined by n, however, it will be necessary 

to select values for the field constants through 
fifth order terms in p to produce a sufficiently accurate 
field. In addition, B may be derived from a potential 
function by the defining equation: 




in p and first order in 5 would have the form: 

"^^31*^ ''’^41^ +f5j^P^+. . . )?] 

at the midp lane is then 

and there is a direct relationship between the field 
parameters b 2 and ^21 



These scalar potential coefficients can be related 
directly to the vector potential coefficients if so 
desired. 

Off-Midplane Field Calculations 

In order to determine the ion trajectories out- 
side the magnetic midplane, it is necessary to ex- 
tend the magnetic potentials to higher order terms in 
?. The vector potential would have the form: 

Ag = A^[ (p+a2p^+a2P^+a^p^+a^p^+. . . ) 

+ (aQ2+a^2P+^22^^'*'®32*^^'^* • • ^ 

^^04"^^l4^*^^24^ +....)§ ] 

where A^ == B^r 
o o o 
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and the B field is determined by 



1 bA0 
=^0 


and 

„ _ 1 bA6 . 1 A0 

®z - + 5;aTpT 

Similarly, for the scalar magnetic potential expan- 
sion to fifth order terms in § we have: 

+(f03+fi3p+f23P^+f^‘^P^+* • • 


.)§ 


■ 33 ' 


+(fo 5 "*'^ 15 P'''* • • 


then 

+3(fo3+fi3p+f23P^+f33p^+- • • 
+ 5 (fQ 3 +fj^ 3 P+...) 5 ^] 

and 

Br = B^[(fj^j^+2f2j^P+3f3j^p^+4f^j^P^+5f3j^P^+. ..)§ 

+(fl3+2f23p+3f33p^+. . .)5^ 

+(f^ 3 +. • 


It can be seen that the above expressions satisfy 
Maxwell's equation for steady state fields, 

V X B = 0 

and by requiring the field to satisfy a second con- 
dition 

V • B = 0 , 
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the off-midplane field coefficients may be deter- 
mined in terms of the midplane coefficients. This 
requirement is 

i)B 

and by equating like powers of p and § these coeffi- 
cients are 

_ (fii+2f2i) 


^03 " ■ 
f = - 

5 

("fll+2f2i+^f32^) 

^13 

6 

f = - 

( f 1 1 - 2 f 2 j^+3 f 3 j^+12 f ) 

^23 

6 

f = - 

( - f 2^ 2^+2 f 2 2^+3 f 3 j^+4 f ^ 2^+2 0 f 3 , 

^33 

6 

f = - 

(-fll+2f2i-12f3i-24f^l) 

^05 

^ 

f = - 

(2fii-4f2i+15f 31-12 f^l-60f 

^15 

60 


These off-midplane magnetic scalar potential field 
coefficients may be related directly to the vector 
potential coefficients by a similar calculation. 
Optical Calculations 

There exist many methods for determining the 
optical properties of a magnetic lens. Point to 
point ray tracing, exact or approximate mathemat\cal 
solutions, and experimental measurements are frequent- 
ly employed. However, the focusing criteria is best 
determined by direct application of fundamental 
theorems. For an optical lens, Fermat's principle 
of least time is applicable, but its mechanical 
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analog, the principle of least action is more 


appropriate. 



is a minimum for any mechanical system where the 
Hamiltonian is conserved or for a particle in a 
central field, 


and 


H = Sp^q^ - L = E (total energy) 





L = T-V 


T = Kinetic energy = %m(f^+r^0^) 

A = Action, V = potential energy. 

Since the Action integral is an extremum, 
tion of the integration path must be zero 
the end points or 
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= 0 . 


any varia- 
even at 


By eliminating the independent variable time, a 

relationship similar to the principle of minimum 

energy is expressed by 
^2 

\ Ep.dq, = K , a constant 
J i a- 1 


and is commonly associated with Maupertuis. Optically, 

this is equivalent to 
2 

^ n*ds = K , a constant 
1 
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where n is the index of refraction. It is easily 
shown that the Hamiltonian formalism and resulting 
equations of motion follow directly from the prin- 
ciple of least action. However, in order to sim- 
plify the analysis, Lagrangian formalism will be 
utilized in the following derivations. 

The relativistic Lagrangian (L) for a moving 
charged particle with velocity V in a static elec- 
tric and magnetic field is given by; 

L = mc^(l-/l-v^/c^) -eU + e/c V*A. 
where U = electric field = 0 for isolated magnetic 
systems 

A = magnetic vector potential 
m * particle mass 
c = speed of light 
e = charge of particle 

For the particular case of the RPI magnetic analyzer, 
using A as previously defined, the momentum components 
are: 


bL 

Pr = -^ = myr 


Pg = myr^e-mYvr (p+a 2 P+- • • ) 

/I 2 / 2n -% 

Y = (1-v /c ) 2 


where by equating the centrifugal and magnetic forces, 
m Ve 


B. 


er. 


also A„ = B^r^ 
o o o 


It is no longer necessary to retain the relativistic 


23 



correction term y since it is negligible and, in the 
least action formalism, it cancels out. However, it 
is essential for relating the particle momentum to 
the magnetic field strength. 

In continuing the investigation of the optical 
properties of a magnetic lens, only the first vector 
potential expansion coefficient (a 2 ) will be retained. 
Also, the basic field inhomogenity (n) is related to 

32 by 

-n = 2 a 2 +l • 

However, the selection of n can be made without re- 
gard to the higher order field expansion terms. The 
final ion trajectory calculations will include terms 
up to fifth order, as derived previously. 

The kinetic energy, T, is given by 
T = %mv^ = %m(t^-h:^0^) 

and 


thus 


~ r 2 2 Ps 2 pgmv^(p+a 2 P^) 

Pr = c® ''o - T 




The initial system conditions are: 
P - 0, r = 

0 = 0 


r = V Sin a 
o 

0 = — Cos a 
o 
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where a is the initial radial angular divergence. 
In terms of p 


Pq = mrQ^(p+l)^0 + er^(P+l)A 


0 


= mr V Cos a 

= mroVo(l 2 + • • •) 

The circular action is then 


I Pgde = mr^v^Cl- |-)§ 


where § = n , for this system. 

The radial action to first order is found by 


?r = ± mv^ |a^[l-p+(2+a2)P^]-2(H-a2)p^}‘ 




c p 

max P max 

. '’r‘*r ' 2 ^ 


Pj.dp 


*•0 o 

where p max may be calculated from the condition 

p = p 

max 

when Pj, = 0 . 


Thus 




max 


max' '• max Z^max" 


or 


= 2(l+a2)p^ 


max 


and the action integral is 


max 


•, % 


2r^ ^ + mv^ |a^[l-p+(2+a2)P^]-2(l+a2)P^| dp 
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factoring out constant terms yields: « 

^ P^max ^ ry [l-p+(2+a^)p‘^] 

2roMv^V2(l+a2) ^ 2(I+a,) ^ } ‘“o 

O ' Z' 

A change of variables of the form 

= f(p)^+p^ where p = X sin |a 


and 


also 


. a^Cl-p+(2+a,)p^] 

TTnsjr — 


2 2 2 
X = ignoring terms higher than p 


thus 

2 




o 


in p space. 

Therefore the radial action is given by 

Action (r) = ^ do 

o 

= mr _ V V 2 ( l+a« ) } 

00 *' 2' 12 maxJ 

The azimuthal action in terms of p is 

max 

Action (0) = mr^v^[l-(l+aQ)p^ ]$ . 

' 00 ' 2' max 

Total action must be independent of or 



for proper focusing. In terms of the field expansion 
coefficient , 

bi - 2a2+l = ^ - 1 
and for the 255° analyzer 

$ = n b^ = -n = -% 
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Thus the basic field inhomogenity of 0.5 should be 
used and the higher order terms may be adjusted to 
maximize radial focusing. By extending this calcu- 
lation to higher order terms, optimum theoretical 
values may be determined for other field expansion 
constants. However, the truncated theoretical 
field parameters do not, in general, provide the 
most ideal focusing conditions. The major intent 
of this study is to determine experimentally, with 
the aid of computed trajectories, optimum field 
expansion parameters through the fifth order. 
Trajectory Calculations 

The theoretical field expansion and equations 
of motion developed in the previous section can now 
be applied to the ion trajectory calculations with 
the aid of a high speed digital computer. Basically, 
the program calculates each ion orbit in a step- 
wise process from its point of origin through several 
hundred steps to its focal point at 255°. Since the 
B field is defined at every point and the force on 
the ion is given by Lorentz as 
F = q(V X B), 

it is only necessary to solve the resulting second 
order differential equations. The numerical methods 
developed by Nystrom^^ were employed for the integra- 
tion in three dimensions. Although the program was 
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designed to specifically calculate alpha particles 
in a 255° system, it can be modified to calculate 
ion trajectories for any sector type magnetic 

analyzer with focal and object point in the field (see Fig. 8). 

In order to obtain the most desirable focus- 
ing properties, the input field parameters listed 
in Table B were varied as indicated. For calcula- 
tional purposes, the field is assumed to extend 
from R = 0 to infinity, however, the calculated 
maximum and minimum radial position of the peripheral 
ion trajectories may be used to determine the re- 
quired field extent for any specified focusing 
critera. It should be noted that several dif- 
ferent combinations of field expansion coefficients 
will produce the focal properties of figures 9 
through 18. In fact, the final values of the higher 
order terms were determined by trial and error 
variance methods. The basic field inhomogenity 
value of 0.5 was chosen so that the source and 
focus points would fall within the magnetic field, 
thus eliminating a detailed azimuthal fringing 
field analysis (see figure 8). However, the radial 
fringing field of this analyzer will be considered 
in order to extend the usable field area as much 
as possible. 

The device's transmission f characteristics is 
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The focal position for 255°, 180°,' 90°, and 60° 
sector magnets using symmetric object-image 
locations. 


Figure 8 
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Radial focusing trajectories of 
the RPI 255® Magnetic Analyzer. 
d= -1.0® and radial angles 
«= +15® to -15® in 3® steps. 




AXIAL DISTANCE (INCHES) 


1.2500 


I .0000 


7500h 


.5OOOU 


T 


• 2500F 


0.0000 


- .2500: 


- . 5000 


- .7500 


1.0000 


t^e^Rpf trajectories of 
fl! ^Ma^netic Analyrer. 

.. 7^0 radial ancles 

®= +15“ to -15“ 


ang les 
3 “ steps. 


-I .2 5001 


_L 


_L 


X 


e46 00 2.7.00 24800 249.00 25000 251 .00 252.00 253,00 254.00 
AZIMUTHAL ANGLE ( DEGREES) 


Figure 16 


3 




AZI 


T 


T 



Radial focusing trajectories of 
the RPI 255® Magnetic Analyzer. 
3= -2.0° and radial angles 
4t = +15° to +15° in 3° steps. 


-I I I I I 1 

250.00 251.00 252.00 253.00 254.00 255.00 

HAL ANGLE (DEGREES) 


Figure 17 



AXIAL DISTANCE (INCHES) 







Table B 


Ion Trajectory Calculation Parameters 


Parameter 

Values applied to 
the 255° system 

Final Design 
Parameters 

Mean Radius (R^) 

25.0 - 30.0 inches 

27.8 inches 

Source 

(R,Z,0) 

Point Source 
30.0,0.0,0.0 

Point Source 
27.8,0.0,0.0 

Particle Mass 
(a - particle) 

4.0 

4.0 

Particle Energy 

5.0 + 0.5 MeV 

5.0 + 0.5 MeV 

Sector Angle 

255° 

255° 

Radial Acceptance 
angle (a) 

-15.0° to + 15.0° 
(3° steps) 

-15.0° to + 15.0' 
(3°steps) 

Axial Acceptance 
Angle (e) 

-2.0° to +2.0° 
(1° steps) 

-2.0° to + 2.0° 
(1° steps) 

Magnetic Field 
Parameters 

-0.5 

-0.5 


0.125 

0.125 


0.0625 

0.0625 


-0.10 to -0.15 

-0.125 

^5 

0.75 to 0.125 

0.095 

Extent of 
Radial Field 

16 to 44 inches 

16 to 44 inches 

Gap width 

3 - 5 inches tapered 

4.0 inches at 


R = 27.8 inches 
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determined by the radial and axial acceptance angles. 
Resolution may be determined by the formula 

Ml = W 
E 'T 

where 

AE = the expected energy deviation of the 
particle due to focusing properties 

E = the particle's energy 

W = focusing image width 

d = the diameter of the magnetic analyzer 
T = a given transmission value. 

The maximum radial acceptance angle is determined 
by the extent of the usable magnetic field. For this 
system, the radial acceptance angle is limited to 15 
degrees. As the figures indicate, the radial focal 
point diverges as radial angular acceptance is in- 
creased. From Fig. 13 it can be seen that for the 
maximum transmission value of 0.37<, of 4tt steradians, the 
corresponding energy resolution is 

AE I _ 1 

E '0.3% T5DU 

By limiting the azimuthal acceptance angle to + 1.0 
degree, the energy resolution is from figure 11, 

AE I _ 1 

T~ '0.157= “ -SUDIT 

while the calculated maximum resolution for midplane 
focusing in Fig. 10 is 

^ I ^ 1 

E 'midplane 12000 
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Thus it can be seen that the selected values 
of field coefficients result in satisfactory optical 
characteristics of the system. The azimuthal cal- 
culations indicate that the image of a point source 
(at maximum acceptance) will be a line. This is as 
expected since, in a double focusing device, opti- 
mum radial focusing criteria will not produce the 
optimum axial focal pattern. This design is based, 
of course, on optimizing the radial focusing properties 
Focal Pattern Computation 

The optimum field expansion coefficients have 
been selected for focusing 5 MeV alpha particles; 
the transmission and focal properties of alpha par- 
ticles with varying energies may be investigated. 

By employing the same method and programs used in 
the 5 MeV case, the focal patterns for alpha par- 
ticles with energies varying from 4.5 to 5.5 MeV in 
0.1 MeV increments can be calculated. The plane 
in which these images formed made a 42° angle with 
the mean radius orbit at an azimuthal angle of 
255°. 

While examining the particle paths from source 
to focal point during this energy sweep calcula- 
tion, it was found that the maximum radial deviations 



from the mean radius were not symmetric. The max- 
imum radial extent of the 5.5 MeV alpha particles 
fell beyond the pole piece outside edge, while the 
4.5 MeV inner trajectories were well within the 
inside pole piece edge as shown in Fig. 19. In 
order to optically center the particle paths, it 
was apparently necessary to reduce the mean radius 
and thus equalize the maximum inside and outside 
radial deviations from the mean radius. By vary- 
ing the mean radius down to a minimum value of 
25.0 inches, an optimum value of 27.8 inches pro- 
duced the centered pattern indicated in Fig. 19. 

The maximum axial extent of all calculated trajec- 
tories must be less than pole piece gap at any 
given radius, otherwise the maximum vertical accep- 
tance angle will be reduced. 

The optical mean radius was therefore changed 
from 30.0 inches to 27.8 inches. A reinvestiga- 
tion of the focal pattern indicated that it remains 
approximately valid for the reduced radius design. 

Only a scale factor on the order of 0.927 is re- 
quired to convert the previously calculated 30 inch 
axial and radial plots over to equivalent 27.8 inch 
plots. Table B lists the final field parameters which 
completely determine the theoretical optical properties 
of the RPI 255° analyzer. It is now necessary to calcu- 
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late a pole piece contour which will reproduce these 
physical properties with a high degree of accuracy. 

Shim Contour Calculation 

It is by no means a simple task to construct 
a magnetic system which will reproduce a given 
theoretical field. Practical problems such as 
material impurities, fringing field distortions, 
magnetic saturation, and construction dimension 
tolerances have always plagued the experimentalist. 

In this work, only the two major design problems 
of saturation and fringing fields are considered, 
with considerable emphasis on the latter. 

If one is willing to give up the outer portions 
of a magnetic analyzer which are distorted by fring- 
ing fields, the pole piece design is greatly simplified. 
Since the magnetic lines of force enter and exit 
the pole pieces normal to their surfaces, the pole 
piece contour may be considered to be an equipoten- 
tial surface. It would only be necessary to shape 
the contour according to an appropriate equipoten- 
tial contour line from the theoretical field in 
order to produce an acceptable magnetic field. If 
this approach were taken for the RPI 255° analyzer, 
the maximum transmission value would be reduced. 

The large yoke system was designed to maintain the 
maximum magnetic induction within the iron components 
to values well below saturation. In addition, sharp 
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points and edges were eliminated from the pole piece 
design since the magnetic flux could easily saturate 
such areas. In particular, the inside and outside 
edges were rounded in order to prevent saturation 
effects from distorting the fringing field regions. 

One of the basic goals of this project is to 
maximize the actual transmission characteristics of 
the device. Therefore, a method must be developed 
to adjust the inside and outside edges of a truncated 
equipotential contour in order to compensate for the 
radial fringing fields. A relaxation type calcula- 
tion appeared well suited for this situation. By 
superimposing a radial cross-section of the analyzer 
over a 0.2 by 0.2 inch grid as depicted in Fig. 20, 
it is possible to solve Laplace's equation numeri- 
cally in order to determine the magnetic field strength 
at each mesh point. The shim contour may then be 
varied until the fringing field effect is eliminated" 
and the required magnetic field shape obtained. 

Boundary conditions were specified as Vm = 10 
at the shim surface and Vm = 0 at the magnetic mid- 
plane. The potential at any given point is calcu- 
lated from the potential at its four nearest neigh- 
boring points. This type of calculation is straight- 
forward, except within the gap of the pole pieces. 

Here, the net spacing is modified by adjusting the 
grid potential values since the actual shim contour 
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does not coincide with normal net spacing. Table C 
lists the infinite field contour coordinates while 
Table D presents the results of comparing the ideal 
(infinite contour extent) field shape and the actual 
potential relaxation calculation for a mean radius 
system of 30 inches. Agreement is relatively good 
except at the radial fringing field, which is as 
expected. 

In general, a negative field error at a spe- 
cific radius indicates that the gap width is too 
large and should be reduced. Conversely, a 
positive error may be reduced by increasing the gap 
width. Therefore, it is a simple matter to elim- 
inate the field error at a single point by experi- 
mentally adjusting the gap width at that point. 
However, there is a strong interaction among 
neighboring points and minimizing the error at all 
points simultaneously is a formidable task. Three 
approaches were taken in order to solve the problem 
and produce an acceptable contour. 

1. Linear Correction Method . The first 
method applied was a process where each shim point 
was linearly corrected relative to the calculated 
error at the point, or 

D'(a) = D(a) + K-E(a) 
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where , 

D(a) = gap width at radius a 

D^(a) = adjusted gap width at radius a 

K = linear adjusting constant 

E(a) = the magnetic field error at radius a. 

After applying this equation to each point across 
the extent of the pole pi“ece, the corrected contour 
shape was coded back into the relaxation program. 

This iteration process was repeated several times 
in an attempt to minimize the calculated errors. 

Upon completing approximately twenty cycles, it 
became apparent that this method was not producing the 
desired results. Two problems were encountered. 

First of all, several shim points were being over- 
corrected during the iteration process, resulting 
in wavelike oscillations (positive to negative) of 
the computed error shape. Reducing the value of K 
damped the oscillations, however, after several 
iterations the shim contour was being "broken up" 
by this process. Specifically, wherever the error 
curve switched from negative to positive, the 
corrective process generated a sharp step or dis- 
continuity in the physical contour. Eventually the 
curve would have assumed a sawtooth form. 

2. Contour Smoothing . Since the first method 
appeared to be adjusting the contour in the right 
direction, it seemed that only a mathematical 
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smoothing process was required. Recalling that the 
original magnetic field is the result of a fifth 
order Taylor expansion, it was felt that a high order 
polynominal fit of the corrected shim data would pro- 
duce results close to the optimum shape. This con- 
clusion is based on the assumption that the optimum 
shape is truly a mathematically well-behaved function. 
However, the application of polynomial fits up to 
twentieth order did not give satisfactory results. It 
appeared that the fitting process so distorted the 
existing shape that the linear correcting method 
could not bend the curve toward the optimum contour. 
Even partial fitting of the inside and outside edges 
to separate contours was unsucessful. 

3. Perturbation Technique . Presently a point 
variation technique is being investigated. Essen- 
tially, this involves moving each shim point by a 
small amount and then calculating the effect at all 
other points. This results in a series of simultaneous 
linear equation of number and order equal to the 
number of shim points. This method should be exact, 
except for multipoint movement interactions. How- 
ever, if the individual point correction is relatively 
small, this effect should be minimal. The major dis- 
advantage of this method is the vast amount of computer 
time required for the calculation. The method is 
mathematically defined by: 
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= matrix of variable coefficients 






error at shim point j without shim 
variation 


= error at shim point j with a varia- 
tion in shim point i equal to A 

A = the distance increment of variation 


X. = the distance variation of point i 
which will reduce the total shim 
error to zero. 

The solution of this system of linear equations is 
applied to the contour as follows: 

D^i) = D(i) + X. 


where is the solution at point i. The corrected 
contour is then coded into the relaxation program. 

The results of this approach are very promising, 
however, a major obstacle to be overcome is the size 
of X^. Mathematically, the value of X^ is unlimited 
and varies over that range. However, the shim con- 
tour is limited by its physical dimensions and can- 
not accommodate point movements over several inches. 

A possible solution to this problem is to limit the 
range of values permitted X^ and then minimize the 
contour error instead of attempting to eliminate it. 

The final shim contour coordinates given in 
Table E were determined by shifting only a few points 
at a time during the iteration process. While it was 
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necessary to continuously readjust the interior 
field to compensate for the few point variations 
performed at the outer edges, the final field 
errors listed in table F for this contour are 
minimal. Further reduction was attempted, but 
calculations indicated the result of reducing the 
outermost field errors more would be an increased 
error in the inside field. Finally, as an over- 
all check of the calculational technique, a 0.1 
inch grid mesh model was developed to compare with 
the results obtained from the 0.2 inch model. 

This fine mesh grid was introduced over the 
inside edge of the larger mesh model. It extended 
from 6 to 26 inches as shown in Fig. 19. Boundary 
values were taken from the large mesh program and 
linear interpolation used to generate the inter- 
mediate point values, except along the shim contour, 
Here, the intermediate values were determined by 
smoothing the contour shape or by taking values di- 
rectly from a 0.1 inch infinite field calculation. 

The results of comparing the two models are given 
in table G. While the results are in good agreement, 
a need for a finer mesh calculation is indicated. 

The minimum mesh size is determined primarily by the 
amount of computer central processing time available 
for the calculation. 

The alpha particle trajectory calculations dis- 
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cussed previously required approximately 0.5 seconds 
of central processing time on a CDC 6600 digital com- 
puter. Approximately 1000 trajectories, were calcu- 
lated and plotted for the optimtun theoretical field 
calculation, or about 500 seconds total time. Each 
infinite field calculation required about 10 seconds 
to reach a convergence value of 1 x 10” . The 0.2 
inch mesh model required 12 minutes per run or about 
360 minutes total central processing time. The one 
0.1 mesh program ran 45 minutes, and calculated only 
a 20 inch segment of the 0.2 by 0.2 inch calculation 
in a 0.1 by Q.l inch mesh grid would require on the 
order of 100 hours of central processing time. Even 
transferring the program to the faster CDC 7600 
machine would not bring the required calculation time 
to less than 10 hours. The choice of 0.2 inch mesh 
model was therefore consistent with the resources 
available to complete this investigation. 

Tables C through G are abbreviated summaries 
of data generated by the computer calculations. They 
include only values for integers, even though data 
is available for 0.2 and 0.1 inch mesh points. 
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TABLE C 


RPI 255° ANALYZER SHIM CONTOUR 
INFINITE EQUIPOTENTIAL CALCULATION 
(30.0 inch radius) 


Radial 
Distance 
( inches ) 

Midplane to 
Shim Surface 
(inches) 

Radial 
Distances 
( inches ) 

Midplane to 
Shim Surface 
( inches ) 

l4.0 

1.563705 

31.0 

2.033611+ 

15.0 

1.5831+08 

32.0 

2.067785 

l 6.0 

I.60I+U01 

33.0 

2.1021+59 

17.0 

1. 626607 

3I+.O 

2.137618 

18.0 

1.61+9965 

35.0 

2.173236 

19.0 

I.67I+I+07 

36.0 

2.209303 

20.0 

1.699878 

37.0 

2.21+5795 

21.0 

1.726318 

38.0 

2.282699 

22.0 

1.753677. 

39.0 

2.319996 

23.0 

1.781903 

1+0.0 

2.357667 

2 h.O 

I.8109I+7 

1+1.0 

2.395686 

25.0 

1.81+07 63 

1+2.0 

2.I+3I+O3O 

26.0 

1.871311 

1+2.8 

2.I+6I+913 

27. 0 

1.90251+6 

l+l+.O 

2.511533 

28.0 

I.93I+I+36 

1+5.0 

2.550593 

29.0 

1.96691+0 

1+6.0 

2.589765 

30.0 

, 2.000000 
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TABLE D 


COMPARISON OF POTENTIAL RELAXATION 
CALCULATION WITH OPTIMUM FIELD SHAPE 
(30.0 inch radius) 


Radial 
Distance 
( inches) 

Percentage 

Difference 

Radial. 

Distance 

(inches) 

Perc entage 
Difference 

16.0 

-14.873687 

31.0 

0.000929 

IT.O 

-3.010685 

32.0 

0.001311 

18.0 

-0.475847 

33.0 

0.001438 

19.0 

-0.074085 

34.0 

0.001513 

20.0 ■ 

-0.012983 

35.0 

0.001565 

21.0 

- 0.003465 

36.0 

0.001368 

22.0 

-0.001887 

37.0 •' 

-0.000080 

23.0 

-0.001569 

38.0 

-0.006712 

2 k . 0 

-0.001452 

39.0 , . \ 

■ -0.033373 

25.0 

-0.001389 

4 o.o 

' -0.134371 

26.0 

-0.001372 

4 l.O ' ' 

' - 0 ; 501079 

27.0 

-0.001355 

42.0 

-1.760605 

28.0 

-0.001273 

43.0 

-5.610616 

29.0 

-0.000920 

44.0 

-14.625993 

30.0 

0.000000 
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TABLE E 


RPI- 255 ° ANALYZER SHIM CONTOUR 
FINAL SHIM CONTOUR CALCULATION 
(27.8 INCH RADIUS) 


Radial 

Distance 

(inches) 

Midplane to 
Shim Surface 
(inches) 

Radial 
Distance 
( inches ) 

16.0 

2.0000 

31.0 

17.0 

1.5550 

32.0 

18.0 

1.6400 

33.0 

19.0 

1.7160 

34.0 

20.0 

1.7420 

35.0 

21.0 

1.7720 

36.0 

22.0 

1.8033 

37.0 

23.0 

1.8353 

38.0 

24. 0 

1.8681 

39.0 

25.0 

1.9018 

4o.o 

26.0 

1.9362 

4 l.O 

27.0 

1.9714 

42.0 

28.0 

2.0071 

43.0 

29.0 

2.0436 

44.0 

30.0 

2,0807 



Midplane to 
Shim Surface 
( inches) 

2.1183 
2.1565 
2.1952 
2.23UU 
2.2jh2 
2.3lh3 
2 . 35^9 
2.3960 

2.k37h 
2.4810 
2.5040 
2. 7350 

2.6000 

2.7500 
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TABLE F 


COMPABISOH OP FIHAL CONTOUR FIELD 
CALCULATION WITH OPTIMUM FIELD SHAPE 


Radial 

Distance 

(inches) 

Percentage 

Difference 

Radial 

Distance 

(inches) 

Percentage 

Difference 

16.0 

-13.961971*' 

31.0 

0.001388 

17.0 

-0.1+1*1270 

32.0 

0.001555 

18.0 

-0.009863 

33.0 

0.001830 

19.0 

-0.006222 

31+.0 

O.OO22I+I+ 

20.0 

-0.001510 

35.0 

0.002718 

21.0 

0.003880 

36.0 

0.003236 

22.0 

0 . 0001+90 

37.0 

0.002795 

23.0 

-0.000922 

38.0 

-0.000361 

21*. 0 

-0.001203 

39.0 

- 0 . 001 + 1 * 61 + 

25.0 

- 0.001223 

1+0.0 

0.005789 

26.0 

-0.001130 

1+1.0 

-0.000689 

27.0 

-0.000782 

1+2.0 

0.021716 

28.0 

0.000207 

1+3.0 

-0.797557 

29.0 

0.000969 

l+l+.O 

- 8.317033 

30.0 

0.001253 
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Table G 

0.1 Inch Mesh Check Calculation 


Radial 

Distance 

(inches) 

Midplane to 
Shim Surface 
(inches) 

Percentage 

Difference 

16.0 

2.000000 

-12.660316 

16.5 

1.500000 

-3.790008 

17.0 

1.555000 

-0.064198 

17.5 

1.870000 

0.344968 

18.0 

1.640000 

0.162830 

18.5 

1.696000 

0.092107 

19.0 

1.716000 

0.039396 

19.5 

1.727000 

0.011949 

20.0 

1.742000 

0.007244 

20.5 

1.756500 

0.007198 

21.0 

1.772000 

0.005498 

21.5 

1.787644 

0.002993 

22.0 

1.803307 

0.001149 

22.5 

1.819198 

0.000199 

23.0 

1.835313 

-0.000108 

23.5 

1.851640 

0.000079 

24.0 

1.868174 

0.000903 

24.5 

1.884913 

0.002996 

25.0 

1.901846 

0.008011 

25.5 

1.918969 

0.019460 

26.0 

1.936275 

0.042748 
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II. THE DETECTION OF SINGLE IONS BY REVERSE BIAS 

P-N JUNCTIONS 

In experiments as early as 1961 principal 
investigator reported on the potential advantages of 
using solid state devices for detecting a beam of 
ions. At that time, however, it was not possible to 
detect single charged particles. Subsequently , 
additional efforts were undertaken to improve circuitry, 
detectors, and minimize corona with the result that 
a one-to-one correspondence has been established be- 
tween the arrival of single ions and output pulses 
that can be registered in a scaler. 

The basic theory of a solid state detector shows 
that a "depletion layer" can be made to function as 
an ionization region for electron-hole pairs. 

The distance to which the depletion layer pene- 
trates into each region of the crystal will be deter- 
mined by a solution of the current continuity equation 
together with Poisson's equation. For the p-type 
region this gives 

(N^ + n - p) 

where N^ is the acceptor impurity concentration, n 
and p are again the local electron and hole densities 
and z is the dielectric constant of the crystal. 

An exact solution of these equations is dif- 
ficult but it is possible to learn a great deal about 
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the junction using a rather drastic approximation, 
namely that throughout the depletion layer n and p 
are small compared with N^. Poisson's equation 
then reduces to 
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Xp being the end of the depletion layer where + n 
- p = 0 (on this model the depletion layer ends 
abruptly) and is the potential at the junction 
interface, x = 0. 

When X = Xp,V - Vq * Vp, the potential dif- 
ference across the p-type region, 
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also the extension into the n-type region is 
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The excess charge in the n-type region is equal and 
opposite to that in the p-type region, for the field 
is confined to the depletion layer in this approxima- 
tion. Hence, 



The depletion layer penetrates the two regions in 
the inverse ratio of their ionized inpurity con- 
centrations, or alternatively in the direct ratio of 


their resistivities measured away from the influence 
of the junction. The depletion layer exists almost 
entirely in the high-resistivity material and for p- 
type material the depletion layer width is approximately 


X 


2 


K V 
2neN^ 



In p-type silicon this gives approximately 

X = 3.2 X 10 ^(pV)^cm 
and for n-type 

X = 5.3 X 10 ^(pV)^cm 

where p is expressed in Q cm and V in volts. 

Basically, the depletion region of a junction 
acts as a solid state ionization chamber with one 
electron-hole pair produced for every 3.6 electron 
volts of energy deposited in the active region of 
the device. However, a major difficulty in the 
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application of p-n junctions to direct heavy ion 
detection is the very short range of ions in solids 
(see Figure 21). As a result, the ion loses most 
of its energy in the dead layer or window of the 
detector and deposits very little in the active 
region. There are other disadvantages of using 
solid state detectors to monitor ion beams directly 
First, for low energy heavy ions, much of the 
energy loss is due to non- ionizing, elastic or 
"nuclear” collisions. Secondly, high ion beam 
currents striking the detector cause degradation 
due to damage of the crystal structure, and un- 
wanted impurity introduction. 

\ 

Since a p-n junction is not efficient in 

detecting ions directly, an alternate method may 

be employed. In this technique, the ions strike 

a dynode ejecting one or more secondary electrons. 

Ion to electron conversion is commonly used for 

ion detection. A conversion electrode offers a 

good combination of high sensitivity and rapid 

response. For example, ion-electron converting 

surfaces, together with an electron multiplier, 

are routinely used in mass spectrometry at count 
6 

rates of 10 counts per second. To increase the 
sensitivity of ion detection systems, a combina- 
tion of ion energy and converter surface prepara- 
tion are chosen to give an average electron yield 
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of from 2 to 10 electrons per incident ion. 

In combination with an ion to electron con- 
version surface or dynode, a single p-n junction 
has a number of attractive features as an ion detector 
for mass spectrometric applications. These are: 

1. compactness 

2. ease of mass fabrication 

3. insensitivity to magnetic fields 

4. high speed (10 to 10 pulses per second) 

5. extremely high efficiency in detection 
of electrons, with a good signal to 
noise ratio 

6. ease of interfacing to data processing 
equipment (as opposed to photographic 
plates) 

7. insensitivity to exposure to air or high 
vacuum 

8. compatability with high vacuum systems 

9. insensitivity to shock or vibration. 

The most important characteristic of p-n junctions 
in this application are their compactness, allowing 
a large number of individual detectors to be used, 
and also the ability of the junctions to operate 
in the magnet’s fringing field, in contrast to 
electron multipliers. 

In present commercial instruments which focus 
a number of masses simultaneously in a plane, the 

*Note: At present a limitation exists with respect 
to output circuitry. 
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most common readout is by means of photographic 
emulsions. It is informative to compare the 
advantages of p-n junctions to photographic 
plates in detail. 

First, p-n junctions can provide a single 
output pulse for every ion that strikes the con- 
verter and produces secondary electrons. For a 
high yield surface, the discriminator can be set 
high enough to give both high efficiency and a 
high signal- to-noise ratio. In contrast, the best 
photographic emulsions require a minimum 10^ ions 
at 10 KeV to produce a detectable line image. 

This fact makes it difficult to detect impurities 
at low levels as in analyzing semiconductor im- 
purities . 

Second, with present p-n junctions, it is 
possible to operate, with suitable electronics, 
at count rates from 1 count per second to 10^ 
counts per second without an appreciable dead 
time correction although the energy resolution 
does suffer. Photographic plates, on the other 
hand, are inherently non-linear and require a 
great amount of precalibration. Multiple ex- 
posures are also generally required. 

Third, while p-n junctions produce pulses which 
may be easily counted or readout for computer analysis, 
a photographic film produces data as an array of lines 



of varying capacity on strips of film or glass. 

These strips must be protected from exposure to 
light and developed, entailing a considerable 
time lag. The strip is then placed in a photo- 
densitometer which transforms the data into 
graphic or digital form. Extensive corrections 
must also be applied for isotopic mass effects 
and film characteristics before qualitative re- 
sults may be obtained. In contrast, a p-n 
junction system in conjunction with a high speed 
computer offers a possibility of real time dis- 
play of the data being collected. (See Fig. 22). 

Finally, a photographic detector requires 
a vacuum interlock to insert and remove the plates; 
also, after insertion the plates must outgas at high 
vacuum, thus increasing the time required for 
analysis. In contrast, the junction array, if 
properly fabricated, should not necessitate opening 
the vacuum chamber for long periods of time. 

For a system employing p-n junctions to be 
efficient and practical, an important requirement 
is the fabrication of a stable, high yield secondary 
electron emitting surface that provides a high 
yield for all ion species. A single detector was 
used to examine a number of high yield surfaces as 
reported in the following section. 
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Figure 22 

Schematic Diagram of 
Real Time Display. 
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The p-n junction system was tested on the 
Cascade Ion Beam Analyzer which has been reported 
in the literature Figure 23 is a schematic 
of the device and Fig. 24 is a photograph. Specific 
details of the manner in which the junctions were 
used in this mass spectrometer are outlined in the 
next section. 
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III. Secondary Electron Yield Measurements for 
Metallic Oxides . 

As an integral part of the p-n junction detec- 
tor effort, it was deemed appropriate to make defini- 
tive studies of several surfaces which might be suf- 
ficiently stable to act as a "converter" of ions to 
electrons. It is obvious that a high secondary 
electron yield is essential if a high signal-to-noise 
ratio is to be achieved in any solid state detector 
scheme. The following paragraphs (substantially 
following the dissertation by John D. Stein) suggest 
that our research effort has indeed been successful, 
and that future work is needed only to complete a 
practical system utilizing a large array of solid 
state detectors. All experiments were carried out 
using the Cascade Ion Beam Analyzer , ^ shown in the 
previous section. 

In the detector region, the mass analyzed 
ions strike the target which is. held at a negative 
potential of 10 to 30 kv (see Fig. 25). The voltage 
feedthrough and plate spacing, however, are designed 
to allow the use of accelerating voltages up to 
50 kv. 

The electron detectors used for the measure- 
ments reported are commercial reverse biased silicon 

2 

surface barrier junctions with areas of 25 to 100 mm . 
The dead layer of these detectors were reported by the 
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Figure 25 



2 

manufacturer to be equivalent to 40 mg/cm of gold 
which is equivalent to approximately 200 angstroms. 

The output of the detector is fed into a low noise, 
charge- sensitive, field effect transistor preamplifier, 
and then into a main amplifier with pulse shaping, 
integration, and differentiation time constants of one 
microsecond. The output is then pulse height analyzed 
and displayed by a 1024 channel pulse height analyzer. 

A block diagram of the detector electronics is shown 
in Fig. 26. With the need for detecting low energy, 
one of the most important considerations is high 
energy resolution and low noise. 

The surfaces investigated were prepared by a 
variety of techniques and were chosen because they 
were presumed to have a high electron yield and pos- 
sible usefulness as an ion electron converter for ion 
beam detection. 

The aluminum oxide surface was prepared by 

the anodic oxidation technique. First, the sample 

was polished to a mirror finish and then washed in 

alcohol and distilled water. The aluminum was anodized 

in a 3 % by weight solution of ammonium citrate. With 

anodizing coefficient of 13.7 angstroms per volt, an 

o 

oxide layer thickness of approximately 800 A was 
produced with an applied voltage of 60 volts. After 
anodization the sample was again rinsed in distilled 
water. 
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BLOCK DIAGRAM OF THE DETECTOR ELECTRONICS 









The berylliiom oxide surface was prepared using 
a 47o weight beryllium copper alloy treated according 
to a process developed by Michijima. The process 
was as follows: The. sample was first heated to 400°C 

in atmosphere for one hour and then slowly allowed to 
cool to room temperature. This treatment results in 
the formation of small inclusions of phase alloy 
which contains 10 to 157o by weight beryllium. Elec- 
tropolishing is then carried out in an 0.5 normal 
phosphoric acid solution. This polishing is carried 
out at room temperature with a current density of 
0.1 ampere-per- square centimeter at an applied voltage 
of from 2 to 3 volts. The sample was then washed in 
deionized water and placed in the vacuum chamber. 

For the titanium oxide, cold rolled 997. pure 
titanium foil was used with the naturally occurring 
oxide layer. The sample was merely washed in alcohol 
and then rinsed in deionized water and then placed in 
the vacuum. 

A typical sample of the spectrum produced by 
the multichannel analyzer is shown in Fig. 27. This 
is a complete piilse height spectrum of the secondary 
electrons produced by a singly charged barium ions, 
with a total kinetic energy of 32.8 KeV, striking an 
aluminum oxide surface. This spectrum was taken 

O 

with a 100 mm Ortec surface barrier detector operated 
at a bias voltage of 105 volts and main amplifier swingle 
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SECONDARY ELECTRON YIELD STATISTICS FOR ALUMINUM 
OXIDE BOMBARDED BY 32.8 KeV BARIUM IONS 



integration and differentiation time constants of 1 
microsecond. The electron accelerating voltage was 
-30 kv and the individual peaks correspond to groups 
containing one to seventeen electrons. From this 
spectrum, a number of important observations can be 
made . 

The data, which is also available in digital 
form on punched tape, allows one to use computer 
plotting routines and computation techniques to 
determine the statistical distribution, as well as 
the mean yield, of the surface as a function of 
ion energy, mass, or charge state. This distribu- 
tion is of particular importance in the design of 
electron and photo-multipliers which up to the 
present has only been determined indirectly by the 
study of pulse height spectra. 

It should also be noted that while the separa- 
tion of individual peaks can be clearly noted, it is 
poorer than would be expected from measurement of 
the system resolution, which was 10 - 12 KeV. Also, 
as the group of numbers increase, the resolution 
of the peaks become smaller. A possible explanation 
for this is the fact that there is a probability 
that an electron might be backscattered out of the 
detector before depositing all of its energy, thus, 
the more electrons scattered out of the detector, 
substantially poorer peak separation at higher groups 
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will result. 

The most probable secondary electron yield was 
measured for two surfaces. It should be emphasized 
that the yield reported corresponds to the peak of 
the electron distribution. The more commonly measured 
yield corresponds to the mean of the distribution 
and has been found in this investigation to be from 
10 to 157o higher in value than the most probable yield. 

The ion beam was produced by evaporating a 
solution containing nitrates and oxides of barium, 
strontium, lithium, and cesium, onto a rhenium fila- 
ment. This filament was then placed in the Mattauch- 
Herzog trajectory where the material was ionized by 
resistance heating of the filament. The ions pro- 
duced were then accelerated and then mass energy 
analyzed by a combination of magnetic and electrical 
fields and then focused on the target. 

1. Secondary Electron Yield for Aluminum Oxide 

The first surface analyzed was aluminum oxide, 
since this surface has been reported to have a high 
secondary yield under ion impact. The results 
obtained for lithium 7, strontium 87, and barium 
138, are shown in Fig. 28. The results for cesium 
are shown in Fig. 29. These results are shown for 
an energy range of approximately 12 KeV to 40 KeV. 

The energy range was limited at the lower end by 
the resolution of the solid state detector for single 
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Incident Ion Energy (KeV) 



electrons . The upper energy limit corresponds to 
the maximum output voltage of the power supplies 
used. 

2 , Secondary Electron Yield for Titanium Oxide and 
Beryllium Oxide 

A second surface studied was titanium oxide. 

The yield for cesium ions is compared with that obtained 

for the aluminum oxide surface and a beryllium oxide 

surface in Fig. 29. The reason that the beryllium 

oxide surface was not investigated further was that the 

yield was so low for cesium that it was not deemed 

suitable for use as an ion to electron converter . 

Theoretical considerations which relate to this 

C19) 

experimental work have been published elsewhere > ' 

However, it is interesting to compare an experimental 
secondary electron yield with a Poissonian distribu- 
tion. In comparing the distribution shown in Fig. 30, 
to a Poissonian distribution with the same mean, nor- 
malized to have the same area as under the experimental 
curve (see Fig. 27) it was noted that the experimental 
curve was much broader than a Poissonian distribution. 
Haussler^^®^ has reported a similar result for 
secondary electron produced by 300 eV electrons on a 

copper-beryllium alloy. A possible explanation is 

(21 ) 

offered in a paper by Baldwin and Friedman. ^ ' In 

this paper it is theorized that a dynode surface may 
consist of small regions with varying secondary 
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A COMPARISON OF THE EXPERIMENTAL SECONDARY YIELD 
STATISTICS WITH A POISSONIAN DISTRIBUTION 




emission coefficients. This inhomogeneous surface 
might lead to an output spectrum consisting of a 
superimposition of Poissonian distributions with 
varying means. An interesting test of this theory 
would be the investigation of a hgihly uniform 
surface such as a vacuum cleaved single crystal. 
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IV. ION SOURCES 


A new technique was conceived and developed 
at this laboratory for thermally ionizing small 
quantities of the chemical elements on single- 
filaments, and for minimizing the production of 

spurious ions that give rise to errors in ^isotopic 

( 22 ) 

ratio measurements. The concept of this new ^ 
single- filament source is to utilize a diffusion 
mechanism in addition to thermal surface ionization, 
in order to obtain a higher ionization efficiency 
and a mass spectrum which is not obscured by impur- 
ities or superimposed molecular spectra. The 
sample, together with its isotopic spike, is en- 
capsulated within an appropriate filament that is 
made of a refractory metal of high work function. 

The filament is then operated at a high temperature, 
which causes the sample atoms to diffuse to the 
exterior surface where a fraction becomes ionized 
by thermal ionization. This approach has proven to 
provide a number of advantages over conventional 
single- filament sources. The results and applica- 
tions of this type of ion source have been published 
by W. G. Myers and F. A. White (Analytical Chemistry 
1861-1864 ( 1969 ). 

Work at this laboratory on ion source develop- 
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merit has also resulted in an additional technique 
for improving thermal ionization efficiencies. 
Although the method employed is primarily electro- 
chemical in nature, it supplements the procedures 
used to obtain the desired characteristics of the 
diffusion- thermal ionization source. By establish- 
ing a thermal ionization filament as the cathode 
( -) in a slightly acidic solution, metallic 
elements in solution can be electro-deposited onto 
its surface. This concept is not original in 
itself, however the choice and order of anode (+) 
electrode material employed does provide desirable 
results . 

The Analytical Mass Spectrometry Section of 
the National Bureau of Standards has established a 
procedure using a platinum - platinum cell to purify 
copper for subsequent mass spectrometric isotopic 
analysis. At this laboratory we have discovered 
that if a rhenium anode is used instead, the anode 
will dissociate quite rapidly when a potential 
(>10 volts) is applied, placing rhenium into the 
solution. Then if a platinum anode is employed, 
both the metallic sample and rhenium ions are 
attracted to and deposited on the filament cathode. 

This technique has two important properties 
relating to the further development of thermal 
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ionization filaments; e.g. (1) most metallic elements 
under analyses can be attracted to the ion source and 
trapped beneath a film of rhenium metal, thereby 
preventing pron 5 )t evaporation of, atoms when the fila- 
ment is heated, and (2) under certain conditions, the 
surface produced by this process can be made very 
porous and rough, with a large surface area. The 
resulting filament properties thus enhance the proba- 
bility that neutral atoms will become ionized upon 
reaching the surface. The method also provides 
means of concentrating specific metallic elements 
from solution, without producing large deposits of 
unwanted material (i.e.. Cl”, N02~>S0^~, etc.) on 
the filament surface. 

This procedure has been applied to the analysis 
of a number of metallic elements at this laboratory. 

It is routinely used on samples which contain sub- 
microgram quantities of the metal to be analyzed. To 
date, improved ion beams of copper, zinc, iron, silver, 
lead, nickel, uranium, potassium, cesium, rubidium 
and lithium have been observed on our mass spectrometers. 
For example, using an Li Cl solution containing a 
total quantity of 0.2 pg Li, the ionization efficiency 
was observed to improve by a factor of 20. 

An alternative thermal-diffusion ionization con- 
cept is presented in Fig. 31 . Here, the primary ion 
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beam impinges upon a thin foil maintained at an 
elevated temperature. Ions are implanted and sub- 
sequently diffuse and become re- ionized, accelerated 
and pass through a second analyzing magnet. The 
scheme would greatly reduce ion beam intensity, but 

a pure atomic spectrum should be obtained. 

An ion source has also been developed which 

has been very successful in supplying alkali metal 
atoms over extended periods of time. This source has 
been used to investigate the energy loss of ions 
traversing thin foils (see Fig. 47). 
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V ENVIRONMENTAL APPLICATIONS OF MASS SPECTROMETRY 
AND THERMAL IONIZATION SOURCES 


We are including the results of two recent 
measurements because of nationwide interests in 
identifying the chemical composition of atmospheric 
particulates. Although such measurements were not 
specified as a technical obligation under this grant, 
we are pleased to present them as they demonstrate 
the very high sensitivity of mass spectrometry, and 
they point up the desirability of such future 
studies on a greatly expanded scale. 

ASSAY OF COPPER IN RAINFALL 

The primary sources of metallic pollutants 
and their rate of introduction into the atmosphere 
is generally known. The natural forces and meteoro- 
logical processes regulating their removal, however, 
are not as well defined, receiving comparatively 
little attention to date. The high sensitivity of 
thermal ionization, mass spectrometry to metal 
analysis was applied at this laboratory to observe 
the dynamic "washout” of copper from the atmosphere. 

Figure 32 shows the copper content of fractionally 
collected rain samples taken during a shower in the 
Troy, New York area, June 1968. The samples were 
collected in plastic containers and spiked with an 
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COPPER CONCENTRATION { >jg/cm ) 
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isotopically enriched nitric acid mixture. The 
analysis was performed on dropsized (10 pi) portions 
of the precipitation whose total copper content 
amounted to only 1 x 10"^^ gm in some cases. The 
results shown in Fig. 32 were found to agree with 
theoretical models describing contaminant wash- 
out mechanisms. Graphical analysis of this data 
indicated the copper concentration in the precipita- 
tion follows the expression: 

C(pg Cu/cm ) = ^ (1 - e ) + — (1-e ) 

where h is in terms of centimeters of rainfall. By 
estimating other parameters of each storm, one can 
calculate the approximate copper content in ambient 
air both prior to and during the period of precipita- 
tion. Such information, along with direct air sample 
analyses allows the determination of atmospheric self- 
cleansing rates, and provides a better understanding 
of the complex mechanisms which establish contaminant 
concentration levels. 

ENVIRONMENTAL MONITORING OF URANIUM 

Also of current interest at this laboratory is 
the analysis of trace elements in urban air environ- 
ments. Although the National Air Surveillance Net- 
work is monitoring urban atmospheres for some 15 odd 
metals, the instrumental methods of analysis presently 
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applied to these samples are not sensitive enough 
to detect those metals which exist in trace concentra- 
tions. One example of the sensitivity and value of 
mass spectrometric analysis to environmental surveillance 
has been demonstrated on air samples taken in June 
1969 from downtown Albany, New York. These samples 
were analyzed for their uranium content by employing 
the isotopic dilution procedure, using thermal ioniza- 
tion source techniques^ ^ developed by the authors. 

The results of this particular analysis are listed 
in Table 1 along with the concentrations of a few 
other metals present in the atmospheres of U. S. 
cities as a means of comparison. 

IMPORTANCE OF ISOTOPIC ASSAYS IN FUTURE STUDIES 

We believe that within the near future, not only 
will the elemental "finger-print” of environmental 
systems be recorded and monitored, but also their 
isotopic composition. In addition to being ideal for 
biological and environmental tracer studies, stable 
isotopic analyses offers a very sensitive means of 
observing human induced changes to natural environ- 
mental systems. A change in the isotopic abundance 
of an element will not only provide a measure of the 
amount of contamination introduced, but possibly its 
source. This technique has already been demonstrated 
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with the element 

Those interested in monitoring the impact of 
the nuclear industry on the environment will find 
the rare earth elements particularly useful for 
such purposes. The rare earth elements are attractive 
here because of their relatively low abundance in nature 
and rather high yield from nuclear fission. Based 
upon a predicted power output of 10^ MW (T) in the year 
1990 from nuclear plants, the yearly production of 
stable fission product elements therefrom would be 
110,000, 58,000, 28 j 000 and 24,000 pounds of neodynium, 
cerium, lanthanum and samarium, respectively. The 
introduction of a small portion of this material 
into the environment would be detectable due both 
to the low levels (i.e., 5 x 10”^^ g Sm/1 in seawater) 
presently existing there and the characteristic 
isotopic abundance difference between natural and 
fission product elements. As an example, the isotopic 
ratio of Ce^^^/Ce^^^ in nature is 7.98, whereas that 
from the fission of uranium 235 is about 0.015. 

Given a body of water containing 12 x 10"^^ g Ce/1 
(seawater), the input of only a small quantity of 
fission product cerium will produce a noticeable 
change in the natural isotopic ratio. In the case of 

3 

a river with a flow rate of 26,000 ft /sec (Hudson 
River at Albany), a nuclear power plant leaking only 
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6 X 10 g/sec into the river would produce an 
isotopic ratio change from 7.98 to 7.5 downstream 
even after complete mixing. 

The need for establishing background concen- 
trations and isotopic analyses of these and other 
ecosystems is apparent. It is here that mass 
spectrometry can provide a unique contribution to 
environmental studies. 
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VI. "STATIC” MASS SPECTROMETRY 


One spectrometer project that vjas quite distinct 
from other technical studies relates to a "static" 
mass spectrometer. This work would appear to have 
special relevance to NASA programs relating to (1) 
the absolute calibration of the partial pressures 
of Xenon and Krypton (2) the use of "static" mass 
spectrometry generally for assaying gases at pressures 
of 10"^^ Torr and (3) the "burn-up" measurement of 
nuclear fuel. With respect to this latter item, it 
will be recalled that for many years "fission track 
analysis" has been useful in determining trace 
amounts of uranium in geological and other samples. 

In this new spectrometric method one relates fissioned 
atoms to partial pressures of the noble gases, Xe 
and Kr, and by an examination of the isotopic spectra. 

The distinction -between a "dynamic" and a "static" 
mass spectrometer system is described briefly. In 
conventional gas analysis, the gaseous sample is 
typically allowed to enter the ionizing region of the 
mass spectrometer source from a suitable gas handling 
system through a variable leak valve. The ionizing 
volume of the electron bombardment source is made as 
small as possible to maximize the ionization of the 
sample compared to the residual gases in the spectrometer. 
Usually the pressure in this region is fairly high 
(10 ^-10 ^ torr) , consistent with space charge limiting 

effects, the sample size, and the vacuum required. 
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in the analyzer section of the spectrometer. For hydrocarbons or 
other gases where dissociation is a problem, the inlet lines and 
source region are usually ten^jerature controlled to provide repro- 
ducible spectra. Due to isotopic fractionation effects through 
the variable leak, isotopic ratios for a given sample may need to 
be corrected in time, as preferential effusion of the lighter 
species depletes it in the sample reservoir. Once the gas enters 
the mass spectrometer, it is pumped away by the system's vacuum 
pumps. There is a continuous flow of gas sample. 

Static mass spectrometry, on the other hand, requires the 
admission of one discrete sample of gas to the spectrometer en- 
closure with no punning on the system thereafter. An ordinary 
electron bombardment source may be used; but it is customary to 
open the ionizing region to the spectrometer, as the sample is now 
one of the residual gases in the spectrometer enclosure. 

The RPI Static Mass Spectrometer is shown in Fig. 33. Figure 34 
is the overall schematic diagram, and Fig. 35 is the ion source for 
producing ionized atoms of the noble gases. The spectrometer 
consists of (1) a 6" radius of curvature magnetic analyzer (2) a 
bakeable analyzer tube with an electron bombardment type source 
(3) a special gas san^le inlet system and (4) an electron multiplier 
detector that operates in the pulse-counting mode. A severe re- 
striction for the multiplier is that it must be able to withstand the 
high temperature bake-out prerequisite to achieving an initial 
analyzer tube pressure of 10”^^ torr. 

The salient features of the burn-up experiment and the sub- 
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Figure 33 "Static" Mass Spectrometer for Analyzing 

Fission Product Gases 
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FIGURE 34 60 SECTOR MAGNET 




DETAILS OF IONIZATION VOLUME 
SHIELD DRAWING OUT FOCUS PLATES 



2000 VOLTS DC 


FIGUIRE 35 ION SOURCE SCHEMATIC 
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sequent assay of fission product gases will now be summarized: 

1. Small samples (~1 ng) of U-235 were deposited on thin 
strips of high purity Al. 

2. These strips were then "folded" so as to stop the 
fission products produced by neutron- induced fission. 

(The mean. range for gross fission products has been 
measured in Al to be about 1.4 x 10”^ cm. (0.00055 in) 

3. The foils were then irradiated in a well known and 
highly "thermal" neutron flux at an nvt of 4 x 10^^ 
n/cm^ or 0.02% U-235 fission burnup. 

4. Upon irradiation, the fission products are assumed to 
be totally retained with the foils. 

5. Subsequent to this irradiation, these foils were 
totally evaporated by a tungsten filament under static 
conditions within the mass spectrometer, thereby 
releasing all fission products including the noble 
gases Xe and Kr. 

6. Ions were then produced by the electron bombardment 
source and isotopic spectra were obtained from the 
electron multiplier/data acquisition system. 

The Xenon and Krypton isotopic ratios were monitored for 
periods up to 8 hours on total amounts of any Xe or Kr isotope of 
approximately 10^^ atoms. Observed ion counting ratio were in the 
range of several thousand ions /sec for the major isotopes. 

The nomograph of Fig. 36 gives the order of magnitude of the 
number of fissions to be expected from an initial sample of U-235 
as a function of integral neutron flux (thermal). Table I gives 
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FIGURE 36 NOMOGRAPH OF BURNUP VERSUS INTEGRATED FLUX FOR U-235 


a computed Krypton and Xenon production for approximately 4.95 x 10^^ 
total fissions, in terms of atoms-isotopic a partial pressure in the 
mass spectrometer. 

Table I 


ISOTOPE 

Kr-83 

Kr-84 

Kr-85 

Kr-86 

Xe-131 

Xe-132 

Xe-134 

Xe-136 



1.440 X 10^° 1.851 X 10"^° 
2.144 X 10^° 2.754 x 10'^° 
3.920 X 10^® 5.039 x 10'^° 
3.197 X 10^® 4.110 X 10"^® 


Figures 37 and 38 show the isotopic spectra produced for one 
sanqjle, promptly after evaporation and after a 7 hour period of 
operation. 

Mr. George Wood of NASA-Langley kindly arranged for the im- 
portant ionization gauge calibration for this work. Krypton and 
xenon calibration gases in pyrex bottles were attached to the 
system through graded glass to metal tubulation. The valves in 
the gas manifold were bakeable, 1/8 inch diaphragm type. 

To pump the system to the ultra high vacuum range, an 80 
liter per second sputter ion pump was selected. To suppress 
hydrocarbon and noble gas backstreaming , a liquid nitrogen trap 
was added above the pump. This also provides an additional pumping 
speed for condensable gases and serves to reduce the fringing 
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FIGURE 37 SAMPLE RC, XE AND KR FREE FROM INTERFERENCE 




field effect from the magnet on the ion pump, by moving the pump 

further from the source - ion beam region. Figure 39 is a block 

diagram of the vacuum system, including the calibration gases. 

Further details relating to this work are described in the 

(25) 

doctoral dissertation of G. Merrill, and the possible implica- 
tions of this technique are being discussed with NASA- Langley . 
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figure 39 BLOCK DIAGRAM OF THE SPECTROMETER VACUUM SYSTEM 











VII NEW DEVICES AND TECHNIQUES 

A. Coincidence Scheme 

A coincidence scheme was devised by the prin- 
cipal investigator (NASA Patent Case No. 7813) in 
which ions are detected. A special feature of the 
invention is that single ions can be detected in a 
radiation environment by causing secondary elec- 
trons to be sensed simultaneously in two electron 
multipliers (or p-n junctions). Figure 40 is a 
schematic of the basic concept. 

B. Microphonic Detector for Neutral Beams 

A new concept has been devised for detect- 
ing neutral particles by sensing the momentum trans- 
fer imparted to a microphone junctioning as a "target 
for a mass spectrometer beam. The inqjlications of 
this concept have been disclosed in a preliminary 
patent disclosure to G. M. Wood, July 14, 1971, (by 
F. White and H. Rasekhi). It is anticipated that 
further work on this project will be undertaken at 
the NASA- Langley Research Center. 

C. New Fabrication Method for Electrostatic Lens 

An improved method for constructing large elec- 
trostatic lenses was devised by J. Stein which re- 
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SCHEME FOR DETECTING SINGLE IONS BY COINCIDENT 
DETECTION OF SECONDARY ELECTRONS 


suited in a publication in the Review of Scientific 


Instruments . This method was used in the fabrica- 
tion of the large electrostatic lenses made for the 
Cascade and Tandem Mass Spectrometer, (see Fig. 41) 
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Figure 41 




VIII INTERACTION OF IONS WITH SOLIDS 


During the period of this grant some effort 
was expended in connection with a NASA program to 
explore an interaction phenomena. For this report 
a summary only of some instrumental aspects of this 
work is presented. Details are available in the 
doctoral dissertation of H. Rasekhi^and the article 
titled "A Tandem Ion Analyzer of Large Radius" by 
Rasekhi and White (Int. Jl. of Mass Spectr. & Ion 
Physics, 8, 277-291 (1972). 

The mass spectrometric facility that was 
designed and constructed comprised two electro- 
static and two magnetic 90° lenses, each with a 
mean radius of curvature of 122 cm. Both the 
electrostatic and magnetic components can be oper- 
ated independently or programmed in tandem, per- 
mitting mass and energy resolution of primary ions 
and secondary species produced in collision pro- 
cesses. A special ion-detector system provides 
"on line" monitoring of the primary ion beam and 
counting of secondary ions and neutral atoms at 
counting rates from 1/sec to lO^/sec. Special 
features of the instrument make it useful in in- 
vestigations of ion energy loss, charge exchange, 
sputtering, and ion implantation, as well as for 
isotopic abundance measurements. 


(27) 
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The configuration of the apparatus is patterned 

after the four-stage mass spectrometer reported by 

(28) 

White and Forman. The size of the instrument, how- 
ever, makes this present apparatus closer to the 
category of isotope separators than to mass spectrom- 
eters. The two electromagnets are represented 
schematically in Fig. 42. Each has a mean radius of 
curvature of 122 cm and an effective angle of 90°. 

The actual physical angular sector is approximately 
86 °; this reduction from 90° is to account for the 
fringing field which extends substantially beyond 
the edge of the pole pieces. 

Figures 43 and 44 are photographs of the 
analyzing system. The electrostatic lenses were 
constructed with a mean radius of curvature of 122 cm, 
thus matching the radii of the electromagnets. The 
fabrication of these lenses was patterned after the 
manner previously described by Stein and White. 

The interelectrode spacing of the lenses is 10 cm. 

With this spacing 100 keV ions can be analyzed 
utilizing a commercially available power supply of 
10 kV. A schematic of the electronics is shown in 
Fig. 45. 

Ions are monitored in a pulse counting 
mode. The output pulses of the electron multipliers 
are amplified in the amplifier discriminator system 
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TANDEM MAGNETIC ANALYZER SYSTEM FOR ATOMIC INTERACTION STUDIES 
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Figure 43 


tMI M*CU*i 
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Figure 44 





BLOCK DIAGRAM OF THE R.P.l. TANDEM MAGNETIC ION ANALYZER SCHEMATIC 

Figure 45 
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which was originally developed by Sawada. Two 
major improvements have been made. The first sub- 
stantially reduced the nonlinearity in the dis- 
criminator circuit. Secondly, an addition of one 
more amplification stage to the preamplijBL er and 
sections improved the input-output characteristics 
of this system. This improved system is now cap- 
able of sensing pulses with amplitudes as low as 
3 mV and effective duration time of 8 nsec. The 
discriminator level can be adjusted between 3 mV 
to 150 mV. Its output is characterized by pulses 
with about 3.5 V in amplitude and 8 nsec in effec- 
tive width. 

It is appropriate to make special mention of 
the simultaneous counting capability for measur- 
ing neutral atoms and positive ions. As indicated 
in Fig. 42 ions becoming neutralized in a collision 
process will proceed through the magnetic field 
and be intercepted by the neutral beam multiplier. 
Positive ions undergoing energy degradation in a 
target (i.e. in a thin- foil) must be focussed onto 
the final detector by reducing the magnitude of the 
magnetic field of magnet No. 2 with respect to that 
of magnet No. 1. There is a distinct advantage 
in monitoring neutrals and ions simultaneously, 
inasmuch as the neutral beam monitor can be used 
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to normalize variations in the primary ion beam. 
Specifically it has been possible to obtain con- 
sistent data at counting rates ranging from 10^ 
atoms per second to 5 atoms per second. 

The data acquisition system is a PDP 8/1 com- 
puter. Output pulses of the discriminator multipliers 
of both neutral and ion detectors are fed into 
individual scalers which are then interfaced to the 
computer. The computer is then programmed to 
furnish ratios, perform other desired mathematical 
operations on the detecting channels, and provide 
a suitable print-out. 

As indicated above the instrument is especially 
useful for investigating atomic collisions and 
molecular dissociation phenomena. As a single ex- 
ample, consider the measurement of atomic stopping 
power of ^Li ions traversing a thin film of nickel. 
Reference is made to the schematic diagram of Fig. 

42 . We should like to determine the energy loss 

incurred by these singly charged ions in traversing 
the foil. Lithium ions are generated in the thermal 

ionization source, accelerated to a kinetic energy 

o 

of 5-30KeV and focused on a very thin (250 A ) foil. 
The actual shape of the energy loss curve is im- 
portant inasmuch as it permits some analysis to be 
made with respect to the nuclear and the electronic 
component of atomic stopping, and transport phenomena 
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29 30 

generally ’ . The nuclear component is essentially 

Rutherford scattering, and the electronic component 
consists of atomic interactions that give rise to 
ionization and excited states of atoms within the 
solid. 

The primary analyzing magnet was used to deter- 
mine the kinetic energy of the lithium ions incident 
upon the foil. The second analyzer was then pro- 
grammed (decreased in field strength) so as to focus 
those ions, which had traversed the foil, onto the 
electron multiplier. Further, lithium ions which 
had undergone charge exchange in the foil and emerged 
as neutral atoms were detected with the neutral beam 
electron multiplier. A ratio was then taken between 
the ion and neutral atom counting rates. This pro- 
cedure eliminated substantially all errors associated 
with variable ion emission from the source, and it 
permitted the acquisition of reliable data at very 
low counting rates. Figure 46 is a single energy 
loss curve of ^Li ions, obtained for a primary ion 
energy of 15.1 Kev. Figure 47 shows some of the 
construction details of the thermal ionization source 
used to provide the beam of primary lithium ions. 
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MATERIAL TO BE 
IONIZED 
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ION SOURCE 






IX. CONCLUSION 


This terminal report indicates a scope of mass 
spectrometric effort that encompasses several im- 
portant areas, e.g., ion optics, the assay of trace 
metals, detectors, ion production, and the inter- 
action of ions with matter. Actually, all of these 
topics are interrelated in the sense that each adds 
to a broader understanding of modern mass spectro- 
metric instrumentation, and collectively they point 
up the ever increasing potential of spectrometry 
in science and engineering. 

It is to be hoped that this effort will be 
resumed at a future date, with a somewhat greater 
degree of coupling than heretofore with the mass 
spectrometric program of the NASA-Langley Research 
Center. There are several specific projects that 
appear to merit an extended research program: 

1. The application of inhomogeneous ion 
optics to an environmental assay field 
instrument . 

2. The development of a microphonic neutral 
beam detector. 

3. Basic studies relating to photon spectra 
produced by ion impact. 

4. An extension of device work on p-n junctions 
and associated circuitry. 
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5. A small spectrometer for monitoring par- 
ticles in the ambient air. 

6. The development of a film detector whereby 
single ions can be registered on an emulsion, 
in a fashion similar to that for recording 
fission tracks. 

7. Basic studies of channeling at low ion 
energies , and 

8. The utilization of "static" gas mass spec- 
trometry in conjunction with ion implanta- 
tion techniques. 
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